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Cross  Field  Transport  of  Electrons: 
Implications  for  the  POLAR  Code, 
Spacecraft  Charging 


1.  INTRODUCTION 

The  emission  of  a  high  power  electron  beam  from  a  vehicle  in  space  can  lead  to  a  positive 
spacecraft  potential  thousands  of  times  greater  than  the  ambient  plasma  temperature.  A  requirement 
of  the  POLAR  code  1  -2  is  that  it  be  able  to  model  this  Interaction  accurately  in  three  dimensions, 
including  the  Earth's  magnetic  field.  Recent  work  at  AFGL  has  been  aimed  at  validating  POLAR  in  this 
limit  by  comparison  with  Magnetron  Theory  as  a  way  of  testing  its  electron  tracking  accuracy.  POLAR 
predictions  of  electron  collection  in  a  magnetized  plasma  will  ultimately  have  to  match  or  be 
reconciled  with  the  theories  of  Rubinstein  and  Laframbolse^  and  Parker  and  Murphy^.  These  studies 
however,  use  either  drift  theory  or  orbit  classification  to  avoid  detailed  trajectory  tracking,  thus  not 
allowing  the  detailed  evaluation  of  the  POLAR  method  that  we  seek. 
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The  planar  magnetron  has  been  studied  by  many  authors;  Page  and  Adams®,  Blrdsall  and 
Bridges®  both  present  analysis  of  greater  detail  than  that  provided  here,  which  is  only  sufficient  to 
properly  identify  the  parameters  that  bear  upon  our  application.  Our  approach  is  then  to  use  a  space 
charge  limited  diode  model,  with  magnetic  field,  to  characterize  the  sheath  about  a  probe  or  spacecraft 
in  the  ionosphere. 

2.  THE  B-SHEATH 

The  analysis  that  we  present  here  is  that  of  a  steady  state  ’Magnetron’  planar  sheath  with  crossed 
electric  E  and  magnetic  fields  B  (Figure  1).  It  is  not  a  complete  magnetized-plasma-probe  theory,  but 
we  feel  that  our  ’B-Sheath’  solution  can  provide  insight  similar  to  that  provided  by  the  Chlld- 
Langmuir^'®  model  for  the  non-magnetized  sheath  by  identifying  the  key  parameters  and 
characteristics  of  the  problem. 

Figure  1  Illustrates  the  situation  where  electrons  are  emitted  from  the  sheath  edge  at  z  =  0.  with 
density  n^,  and  velocity  z^.  We  take  the  transverse  velocity  to  be  x^,  =  0  at  the  sheath  edge,  as  well  as  the 
potential.  V.  The  electron  charge  and  mass  are  -e  and  m.  B  is  in  the  p  direction  and  E  Is  in  the  ^ 


direction. 

Our  equations  are  thus; 

Force. 

mx  =  eBz  ( 1) 

m2  =  eE  -  eBx  (2) 

Energy, 

jm(x2  +  z^)  -  eV  =  ymZo  (3) 

The  connection  between  a  diode  model  and  a  plasma  sheath  Is  made  by  recognizing  that  the  flux  is 
provided  by  some  net  drift,  v^,  so  for  continuity  we  have, 

n(z)z(z)  =  n^vd,  (4) 

and  Poisson’s  Equation. 

d^V/dz^  =  en(z)  =  (en^Zo/eoZ)  (5) 


5.  Page.  L.  and  Adams.  N.I.,  Jr.  (1946)  Space  charge  in  plane  magnetron,  Phys.  Rev.  69:9.  10. 
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Figure  1.  Diagram  Illustrating  a  crossed  field  diode  or  plasma  sheath.  The  sheath  edge 
Is  modeled  as  a  cathode  where  both  the  electric  field  and  potential  are  taken  to  be  zero. 
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where  Eq.  (4)  has  been  employed  and  we  have  defined  =  v^/z,,  as  the  electron  mach  velocity  by 
identifying  z^  as  the  electron  thermal  speed. 

Integrate  Eq.  (1).  set  Xq  =  0.  and  substitute  into  Eq.  (3)  to  get 

coj^^z^  +  -  2eV/m  =  Zq  (6) 

where  coc  =  eB/m  is  the  electron  cyclotron  frequency.  Differentiating  twice  by  z  (using  d/dz  =  z~’d/dt) 
gives, 

+  z/z  -  (e/mld^V/dz^  =  0  (7) 

Combining  Eqs.  (5)  and  (7)  we  have 

Z  +  (Oc^Z  -  Wp^MeZo  =  0  (8) 

where  Wp  =  yn^eVeom  .  is  the  plasma  frequency. 

The  solution  to  Eq.  (8)  is 

Z  =  ZoMc'Ml-M^.COp^/COc^lSin  +  (McC0p2/(i)c^)Zot 

Define  q  =  cop^/coc^  as  the  characteristic  parameter  with  qj  =  Meq  as  an  effective  q.  Identify  a  =  mzo/eB 
as  the  ambient  gyro-radius,  and  put  the  solution  in  the  normalized  form, 

Z  =  z/a  =  (1-qa)  sinwct  +  qa^ct 

with. 


Z  =  z/COc  =  (l-Qd)  COSCOct  +  Qcl 

2 

Define  a  normalized  potential  4>  =  V/(mzo/2e),  and  return  to  Eq.  (6)  to  get 

<t>  =  (Z2  +  z2-  1)  (9) 

<t>  could  be  obtained  as  a  function  of  t,  but  with  an  unnecessary  Increase  In  algebraic  complexity,  since 
we  can  still  obtain  <I>  as  a  function  of  Z  by  obtaining  Z  and  Z  from  t.  Note  also  that 

q  =  (0p2/a)c^  =  nm/(eoB^)  =  (10) 

where  =  Ve^kT/n^  is  the  Debye  length. 

Figure  2  shows  the  results  of  our  calculation  for  selected  values  of  q.  The  calculations  were 
halted  when  Z  was  observed  to  change  sign,  indicating  that  the  electron  was  turning  around  to 
complete  its  cycloid.  At  this  point  it  is  helpful  to  make  a  distinction  between  a  transmitting  sheath 
where  the  electron  flux  crosses  the  sheath  completely,  and  a  non-transmitting  sheath  where  the 
electron  flux  penetrates  a  limited  distance  (and  corresponding  potential)  and  returns.  The  key  result 
of  our  analysis  is  that  for  qj  >  0.5,  the  sheath  is  transmitting  and  the  total  sheath  potential  drop  is  an 
unlimited  function  of  the  sheath  thickness.  Figure  3  presents  the  same  results  as  Figure  2,  but  with  the 
coordinate  Z.  renormalized  to  the  Debye  length.  From  Figure  3,  one  can  easily  notice  that  the 
magnetic  field  does  little  to  affect  the  sheath  structure  provided  it  is  not  so  strong  as  to  cut  off  the  flow 
(cici  <  0.5).  Figure  4,  taken  from  Page  and  Adams  (1958)  illustrates  the  same  result,  where  the  current 
collected  across  a  diode  is  little  affected  by  B  until  cut  oIT  at  B  >  B^.,  where  Be  is  found  from  Eq.  ( 10)  by 


4 


Figure  3.  The  same  data  as  Figure  2  but  with  the  sheath  thickness  normalized  to  the 
Debye  length.  Values  of  qj  are  Indicated  along  with  the  terminal  points  for  q^j  <  0.5. 
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Figure  4.  A  plot  of  current  vs  magnetic  field  strength  where  the  current  is  normalized  to 
the  Child-Langmuir  space-charge  limited  current  for  B  =  0,  and  B  is  normalized  to  a 
critical  value  corresponding  to  q  =  0.5  (from  Page  and  Adams.  1958). 
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(jj  =  0.5.  For  qj  <  0.5  the  sheath  is  transmitting  only  If  tlie  potential  and  thickness  are  less 
than  limits  obtained  by  setting  Z  =  0,  wliich  gives  oy.t  =  arccos(J-qd'')"'.  For  qj  <  0.5,  this  relation  has 
rcjots  of  it/2  <  wyt  <  k:  for  qj  >  0.5,  Z  >  0  always.  For  a  given  qa  <  0.5,  the  limiting  e),  !  gives  the  limiting 
ih.  iind  Z,  which  arc  ah\'a>'S  small,  i.e.,  Z,„;ix  <  n/2.  <J>  <  Iji/2}^-1  =  1.56. 

Consider  a  non-transmitting  sheath  with  <£>  >  >  <l>n,ax  (g  <  0.5).  For  a  strictly  ID  problem,  these 
conditions  require  a/i  evacuated  gap  large  enough  to  lower  the  capacitance  of  the  gap  to  where  the 
liiniied  charge  of  the  non-transmitting  sheath  equals  that  of  the  surface.  This  is  of  course  a  contlitio'i 
of  niugnrlic  insuhilion  which  can  occur  for  some  gcomelrical  arrangements.  After  consideration  of 
this  possibility,  and  modeling  in  3D  with  POLAli,  which  is  discussed  in  the  following  section,  we  are 
led  to  the  conclusion  that  a  magnetical  insulating  sheath  would  be  unlikely  in  any  space  plasma- 
object  interactions,  and  that  the  same  high  B  conditions  that  produce  a  non-transmitting  sheath  will, 
(hie  to  niinintal  increases  in  cycloid  radii,  allow  charge  to  transport  into  the  olfierwise  evacuated  gap 
iK^m  ends  any  distance  away.  Staled  otherwise,  any  potentially  non-transmitting  sheath  must  be 
considered  first  in  three  dimensions  before  any  simplifications  are  attempted.  It  should  be  noted  that 
in  the  low-earih-orbit  ionosphere,  q  >  1  typically  for  ions,  and  q  <  1  typically  for  electrons. 

A  possible  weakness  in  (his  model  is  the  assumption  (hat  particles  enter  the  sheath  norma!  to 
the  sheath  edge.  Elsewhere.  (Dubs  &  Cooke)^  we  have  reworked  our  analysis  for  enlr>'  at  an  angle  a 
with  t  he  .X  axis  (siieatfi  edgej.  We  find  (hat  (he  minimum  value  of  q  for  which  the  particles  transit  the 
sheath  is  cp,,  =  1  /(2sinn),  I  lere.  we  have  assumed  a  =  90'.  Tliis  model  (and  POIAR)  should  be  augmented 
!)y  a  [jrcsheath  model  tli;U  accounts  for  (he  variation  of  «  with  the  gj-Tation  phase  of  particles  entering 
the  sheath. 

3.  POLAR  RESULTS 

To  lest  POIAR  against  our  B-shcalh  model,  we  constructed  a  flat  octagonal  disk,  pseudo  radius  of 

з. 0  m(;ter.  moving  so  that  the  plasma  impinges  norma!  to  its  face  Vv-iih  B  parallel  to  its  face.  The  disk 
i.s  at  a  positive  potential  r)f  43  volts,  and  moves  at  a  speed  of  mach  8  with  respect  to  oxygen  ions,  and 
nun  h  O.O-ty  will!  respect  to  electrons,  dlie  plasma  parameters  are  0.2  eV  temperature  for  ions  and 

cli  (  Irons,  an  ambient  density  of  2.0  x  lOVem'b  and  a  Debye  length  of  0,74  cm. 

nic  I’DIaR  method  of  elec  tron  tracking  begins  at  a  sheath  edge,  located  as  an  eriui-polential, 

и. si.'.iily  0.09  kT.  External  to  this  siirfac  c  the  electron  distribution  is  presumed  to  be  a  maxwellian 
c  oiistrained  to  flow  only  along  B.  Tlie  dot  product  of  the  local  sheath  surface  normal  with  B  thus 
(I'lcrminc.s  the  flux  and  entry  velocily  of  a  .super-electron  which  is  then  tracked  through  the  .sheath  to 
(ii'icrmine  both  space-charge  density  and  surface  currents.  The  tracking  algorithm  has  an  eHiciency 
nioi  ic.iti  d  .switch  between  direct  integration  of  the  I^rcntn  cciuation  of  motion,  and  gukk'd  renter 

ii  .ijri  Ifiiy,  Tlu-  switcti  to  guiding  center  tracking  is  made  if  it  appears  that  the  radius  of  cui-valure  (not 
.imbii  n(,  but  the  aei  rieralinn  affected  cureaturel  will  be  less  than  the  dimensions  of  the  upcoming 
volume  element.  We  rlemonslrale  that  this  causes  a  grid  dependence  of  the  results  when  the  resolution 


0.  Dubs.  C.W.  and  Cofike,  D.L.  ( 1987)  Particle  Trajectories  and  Potentials  in  a  Plane  Sheath  Moving 
in  a  Magnetoplasma.  AFGL-TH-87-0225.  AD  A196228 
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is  too  coarse,  but  witli  surficieiitly  fine  ^ridding,  gives  remarkably  good  agreement  with  cur  n-shcath 
predictions. 

Model  I  Ls  a  I’niAK-disk  model  with  n  =  0.2  Gauss,  q  =  56.  and  qj  =  2.6.  Figure  5  presents  2-D  cut 
in  the  plane  of  the  flow  vector  (+Z,  from  the  left)  and  B  (+Y.  out  of  the  page.  thus,  along  the  leading  face) 
any  K  x  B  drift  will  be  +X.  Figure  5  shows  the  equl-  potential  contours,  while  Figure  6  shows  the 
sheath  surfaces  and  particle  tracks  with  tfie  contours  removed  for  clarity.  A  wake  structure 
characterized  by  weak  negative  potentials  caii  be  observed  in  the  figures.  Note  that  the  electron 
emitting  sheath  edge  (-  symbol)  is  smooth,  and  most  of  the  electron  tracks  traverse  the  sheath  to  the 
fat  e.  Ttiis  is  tlie  anticipated  behavior  as  both  ()  and  q,i  are  greater  than  0.5.  Figure  7  shows  the  !’()1A]< 
solution  for  B  =  0.  which  can  be  seen  to  be  (lualltatively  similar  to  model  1. 

For  iiKxiel  2.  Figure  8,  B  is  iiu  icased  to  0.3  Gauss  with  q  =  22.0.  and  (la  =1,1.  'I'he  sheath  still  has 
tlie  nominal  no-B  appearance,  but  many  electrons  are  starting  to  drift  along  the  sheath  edge. 

In  model  3.  Figure  0.  B  is  inerea.sed  to  0.5  Gauss,  with  q  =  8.2  and  cia  =  0.38.  I^adically  different 
electron  behavior  is  observed,  with  most  ram  electrons  drifting  parallel  to  the  sheath  edge.  gi\'ing  rise 
to  strong  spatial  perturbations  of  the  sheath  edge.  It  should  be  noticed  that  there  is  still  significtmt 
space  charge  in  the  sheath,  but  that  it  is  now  coming  from  the  edges  of  the  disk. 

Figure  10  Is  a  reiday  of  model  1.  but  witli  the  grid  resolution  reduced  by  half.  'Hie  tendency 
towards  liremature  guiditig  cetiter  treatment  can  be  observed  in  the  increased  liumber  of  drifting 
trajectories  as  if  (la  had  been  ri-duced.  In  this  case,  the  ram  side  sheath  is  resolved  with  about  three 
mesh  intervals.  whi<  h  would  be  the  limit  of  POlAR’s  accuracy  with  or  without  a  magnetic  field  due  to 
the  non-linear  nature  of  ;t  sheath. 

4.  DISCUSSION  AND  CONCLUSIONS 

Linson**^  has  previously  observed  the  importance  of  q  in  determining  electron  sheiith  [iliysics. 
.ititl  the  IciKk-ncy  forcro.ssed  field  (F  and  B)  beams  to  be  unstable  for  ()  -  0.5.  POl-AI^  is  not  a  time 
deiicndent  simulation  code,  but  the  standing  waves  on  the  surface  of  our  lowest  (]  sheath  might  be 
(  onslru<'d  to  be  an  artifai  tual  analog  of  the  electrostatic  instability  invoked  by  hinson  and  studied  by 
Biinemann  and  Linson*  *,  Our  study  of  the  B-shealh  leads  us  to  disagree  with  hinson  on  the  i)oint  that 
lor  (|  j  0.5,  turbulence  is  not  needed  for  the  transport  of  electrons  across  B,  However,  instabilities  at 
the  sheath  edge  m.iv  still  jjlay  an  important  role  in  increasing  the  electron  flux  through  the  sheath 
edge  .md  nariow  the  gaj)  between  the  pla.sma  (|  and  the  effective  q,i  for  the  sheath. 

We  may  now  also  lie  together  our  previous  conclusion  about  the  inescapability  of  end  effects  for 
I  111-  long  evlinder'  problem,  and  (he  problem  (hat  arises  with  qq  when  an  object  is  stationary,  and  the 
ambirni  gyro  r.uliiis  is  small  compared  to  the  cylinder  radius.  In  this  case,  electrons  find  their  way  to 
I  he  slu  alh  edge  only  along  B.  so  there  will  be  no  flux  through  the  sheath  if  its  edge  is  exactly  parallel  to 
1!  riiis  is  no  [jroblem  if  B  is  so  strong  (<j  ■  <  I)  that  charge  trans|)orting  along  B  within  the  sheath  (from 


K)  hinson,  R  H.  ( li''-'))  Current -voltage  characteristics  of  an  electron-emitting  satellite  in  the 
Ionosphere,./,  (,<nf//ii/s.  Rrs,.  74:9. 
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Figure  5.  POLAR  2D  slice  plot  through  the  middle  of  the  3D  grid  and  disk.  Plasma  flows  from  the 
left.  Electron  trajectories  are  started  (symbol  -)  just  inside  the  sheath  edge  (symbol  P).  The  zero 
potential  contour  is  indicated  by  a  0,  and  T  marks  the  negative  potential  sheath  edge  where  ion 
trajectories  could  be,  but  were  not  initiated.  Coordinate  scale  Is  In  grid  units  of  0. 18  meters. 
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MODEL  2  q  -  2Z0  .  qd  -  11  .  B  -  0.3  G  ® 

Mi  -  8.0  .  Me  -  0.05  .  Vo  -  43.0  V  .  R  -  1.8  m 
No  -  20x10«  cc”"*  .  kT  -  02  eV.  Xd  -  0.74  cm 


Figure  8.  POLAR  2D  slice  plot  through  the  middle  of  the  3D  grid  and  disk.  Plasma  flows  from  the 
left.  Electro*  trajectories  are  started  (symbol  ^  Just  Inside  the  sheath  edge.  The  symbol  +  marks 
where  the  ions  are  started  Inside  the  negative  potential  sheath  edge.  Coordinate  scale  is  in  grid 
units  of  0. 18  meters. 
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Figure  9.  POLAR  2D  slice  plot  through  the  middle  of  the  3D  grid  and  disk.  Plasma  flows  from  the 
left.  Electron  trajectories  are  started  (symbol  -)  Just  Inside  the  sheath  edge.  The  symbol  +  marks 
where  Ions  are  started  inside  the  negative  potential  sheath  edge.  Coordinate  scale  Is  In  grid  units 
of  0. 18  meters. 
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Figure  10.  Rerun  of  Model  1,  but  with  Grid  Resolution  Halved  to  0.36  Meter/Grid  Unit. 
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the  ends)  provides  charge  for  the  shielding.  For  q  <  0.5,  however,  this  end-effect  charge  will  still 
undergo  an  enlargement  of  the  individual  cycloids,  causing  some  of  this  charge  to  Impact  the  cylinder 
surface.  This  would  lead  to  an  enlaigement  of  the  sheath  and  a  sheath  edge  making  an  angle  with  B. 
thus  allowing  a  flux  to  the  sheath  edge  and  a  non-zero  qa. 

We  conclude  by  commenting  on  the  value  of  modeling.  In  this  study,  our  end  effect  conclusion 
was  greatly  influenced  by  its  unavoidable  presence  in  the  POLAR  models. 
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